Defect Junctions and Domain Wall Dynamics 
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We study a number of domain wall forming models where various types of defect junctions can 
exist. These illustrate some of the mechanisms that will determine the evolution of defect networks 
with junctions. Understanding these mechanisms is vital for a proper assessment of a number 
of cosmological scenarios: we will focus on the issue of whether or not cosmological frustrated 
domain wall networks can exist at all, but our results are also relevant for the dynamics of cosmic 
(super)strings, where junctions are expected to be ubiquitous. We also define and discuss the 
properties that would make up the ideal model in terms of hypothetical frustrated wall networks, and 
provide an explicit construction for such a model. We carry out a number of numerical simulations 
of the evolution of these networks, analyze and contrast their results, and discuss their implications 
for our no-frustration conjecture. 



I. INTRODUCTION 



If our current understanding of particle physics and 
unification scenarios is correct, defect networks must nec- 
essarily have formed at phase transitions in the early 
universe [1]. The type of defect that forms and its spe- 
cific properties will depend on the particular details of 
each symmetry breaking, so a wide range of possibil- 
ities exist, with correspondingly different cosmological 
consequences — see [2] for a thorough discussion of these 
possibilities. One example is what happens when two ex- 
tended defects (cosmic strings or domain walls) interact 
with each other. In the simplest models they will inter- 
commute (or exchange partners); for cosmic strings, for 
example, this mechanism is crucial for the existence of the 
well-known linear scaling solution (commonly referred to 
as scaling). 

It is possible, however, that topological reasons will 
prevent intercommuting, and that the result of interac- 
tions is the formation of junctions (or bridges) between 
the original defects. An interesting question is then what, 
if any, is the effect of this mechanism on the dynamics 
of the defect networks. In particular, one needs to un- 
derstand the detailed behavior of the defect junctions, 
and the role they play in the evolution and cosmolog- 
ical consequences of the network. These models have 
been much less studied [3-9], though this is purely due 
to the inherent difficulty of doing so, and not due to any 
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lack of physical motivation for doing it. In fact, quite 
the opposite is true: recent work has shown [10-12] that 
cosmic (super)strings will necessarily form at the end of 
brane inflation, and it is believed (though not yet un- 
ambiguously proved) that such networks will have junc- 
tions. Understanding their role should therefore be a key 
step in studying the evolution of these networks, and in 
particular for determining whether or not a linear scaling 
solution is also an evolutionary attractor in this case: de- 
spite some recent supporting evidence [7, 8, 13, 14], this 
remains an open question. Even though understanding 
the evolution of cosmic (super)string networks is not the 
main motivation for the present work, it is clear that 
some of our results provide strong hints for what should 
happen there. 

The present work is the second in a series whose driving 
goal is to understand the cosmological evolution of do- 
main wall networks [15-18] in sufficient detail to be able 
to answer the question of whether or not a frustrated 
defect network [19] can exist at all. In a previous work 
[20] (hereafter referred to as Paper I) we have shown that 
the answer is trivially no for the case of cosmic strings, 
and provided arguments leading to the conjecture that 
the answer should also be no for domain walls. We have 
also established the simple but crucial (and yet previ- 
ously neglected) result that in models with more than 
two vacua, if all the domain walls connecting the various 
vacua have equal energies then only Y-type stable junc- 
tions will form. Here we focus on the role of junctions, 
and we will be studying several models where different 
kinds of junctions can appear. 

Specifically, we consider a number of realizations of the 
two classes of models introduced in [5, 21] (the former has 
also been subsequently considered in [22]). The similar- 
ities and differences between them will allow us to point 
out the key mechanisms at play, and further characterize 



the differences between models with stable Y-type junc- 
tions, models with stable X-type junctions, and models 
where both types can co-exist. The details of the nu- 
merical simulations that we shall present are analogous 
to those in [15-18]. We then go beyond these particular 
models, and discuss the properties that would make up 
the ideal model, at least from the point of view of ob- 
taining hypothetical frustrated networks. We provide a 
general and explicit construction of this ideal model, and 
again numerically study some examples. The pentahe- 
dral symmetry model discussed in [23] turns out to have 
some of these properties of this ideal model. In the in- 
terest of clarity, we shall start by discussing the simplest 
models, and gradually move on to more elaborate ones. 
In particular, we start by considering two-field models, 
then move on to three-field models, and finally discuss 
the ideal model and its relation with other existing mod- 
els. Also in the interest of simplicity, our simulations 
are done in two spatial dimensions, which is sufficient 
to illustrate the key processes. The cosmologically more 
interesting (but also more complicated) case of three spa- 
tial dimensions will be discussed in subsequent work. 



II. TWO-FIELD MODELS 

We will start by considering models described by two 
real scalar fields, which is the minimum configuration 
required in order to form networks with junctions. We 
will discuss in some detail the Bazcia-Brito-Losano model 
(henceforth referred to as the BBL model) introduced in 
[21]. Its original purpose was as a driver for a dynamical 
mechanism leading to 3 + 1 space-time dimensions, but 
from our point of view this particular motivation is not 
essential and we can simply take it as a toy model, which 
turns out to contain a simple yet rich and illuminating set 
of possibilities. The model has the following Lagrangian 
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are real scalar fields and the potential has 
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(2) 
where r and e are two real parameters. This potential has 
minima at the vertices of a square in the plane (<pi, fa), 
the orientation of which depends on the value of the per- 
turbation parameter e that can be assumed to vary in the 
range — 2 < er 2 < 1. There are a total of six indepen- 
dent topological sectors connecting the different minima. 
In the range —1/2 < er 2 < 1 the minima are 
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This model allows for Y-type and X-type junctions de- 
pending on the value of e, which controls the tension of 
the walls connecting each pair of vacua. Note that there 
are two classes of walls which we will denote simply by 
edges and diagonals. In the former the wall joins two 
neighboring minima in field space, and there are four such 
walls. In the latter the wall joins two opposite minima 
in field space, and there are two such walls. In Ref. [21] 
the authors studied only the static behavior of the model 
in the limit for very small e and also made the choice 
r = y3/2. In this case, the minima of potential are in 
the branch given by Eq. (3), that is 4> 2 — (2/3 — e)~ 1//2 . 
By integrating the T® component of the stress-energy 
tensor across the wall between each pair of vacua one 
obtains [21] 
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where a e and ad respectively denote the tension of the 
walls of edge and diagonal type. Depending on whether 
the diagonal walls have a tension smaller or larger than 
twice that of the edge ones, the formation of a Y-type 
or X-type junction will be favored on energetic grounds 
(see Fig. 1). In other words, in the case ad < 2cr e (which 
corresponds to e > 0) we expect to have only Y-type 
stable junctions, while for o~d > 2c e (which occurs for 
e < 0) we expect that only X-type stable junctions will 
be formed. 

Of course, this energetic argument is valid generically, 
not just in the BBL limit e — ► 0. We have analyzed 
the relation between tension of edges, a ei and diagonals, 
ad, for the whole range —2 



< er < 1. We found two 
values of er 2 for which ad = 2er e : in addition to the 
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BBL case, er — 0, there is also 

leads us to expect that for er 2 > and er 2 < — 1 Y- 
type junctions are favored while for — 1 < er 2 < X-type 
junctions are preferred. These solutions are illustrated 
in Fig. 2. In passing, we note that the parameter e 
determines not only the ratio of the energies in the two 
sectors, but also influences, among other things, how fast 
the unstable junctions will decay into the stable ones (this 
is an interesting and relevant question in itself, which is 
left for detailed study in the future). 

These expectations can be confirmed numerically. This 
model is very easy to simulate; in what follows we will 
present the outcome of a few 256 2 simulations. Although 
these are relatively small for today's standards, they are 
more than sufficient to confirm our analysis. For easier 
comparison with BBL, we will also take r = a/3/2. In 
this case, according to the above solution, for —1/3 < e < 
2/3 the vacua are (± (2/3 - e) ~ x / 2 , ± (2/3 -e)- 1 / 2 ), while 
for -4/3 < e < -1/3 they arc (0, ±(2/3 + e/2)- 1 / 2 ) and 
(±(2/3 + e/2)~ 1 / 2 ,0). The outcome of two such sim- 
ulations is illustrated by the four snapshots displayed 




FIG. 1: Depending on whether the thick wall has a tension 
smaller/larger than twice that of the lower tension ones the 
formation of a Y-type/X-type junction will be favored on en- 
ergetic grounds. 
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FIG. 2: Types of junctions (bottom) and configuration of the 
minima (top), as a function of the parameter er 2 in the two- 
field BBL model [211. 



in Figs. 3 and 4, where we have respectively taken 
e = ±0.2. We caution the reader that a pair of Y-type 
junctions close together might be misidentified as an X- 
type junction by a careless glance at a low-resolution 
snapshot. Notice that indeed in the former case only 
Y-type junctions survive, and any X-type junction that 
happens to form is unstable and will decay quickly, while 
in the latter case it is the X-type junctions that survive. 

Again, our treatment is not limited to the asymptotic 
case e — > 0, but it is completely general. Our simulations 
confirm that for any chosen value of e > (that is in the 
range < e < 2/3, given our choice of r), the behavior 
of the system is similar to that in Fig. 3: the formation 
of Y-type junctions is always preferred. On the other 
hand, we confirm our expectation that two distinct cases 
appear for e < 0. For the interval —2/3 < e < 0, only 
X-type junctions are formed, as can be seen in Figs. 4 
and 5 for e = —0.2 and e = —0.4 respectively. However, 
for e < —2/3 X-type junctions are no longer stable, and 
instead Y-type junctions form, as illustrated in Fig. 6 for 
the case e = —0.8. 

As further confirmation, it is also interesting to ana- 
lyze the behavior of the phase space formed by the two 
scalar fields, (pi and fa- Figures 7-10 show the field con- 
figurations corresponding to the snapshots of Figs. 3-6 
respectively. It is quite clear that in the cases e > and 
e < —2/3 the field values tend to concentrate along the 
solutions of both the edges and the diagonal sectors of the 
square — notice that the orientation of the square is differ- 
ent in the two cases, corresponding to the solutions given 
by Eq. (3) and (4) respectively. On the other hand, in 
the —2/3 < e < range the diagonal sector of the square 
is a repeller (which gradually becomes depopulated) and 



FIG. 3: The matter-era evolution of the domain wall net- 
work model of [21] with potential of Eq. (2) and a parameter 
choice r = y/3/2 and e = 0.2. The simulation starts with 
random initial condi tion s insid e the square whose vertices are 
the vacua (i-^/15/7, i-^/15/7). Note that only stable F-type 
junctions survive. From left to right and top to bottom, the 
horizon is approximately 1/16, 1/8, 1/4 and 1/2 of the box 
size respectively. 




FIG. 4: Same as Fig. 3, except that now e = —0.2. 
Here the simulation starts with initial conditions cho- 
sen random ly i n the square whose vertices are the vacua 
(i-^/15/13, ±y 15/13). Note that in this case only stable X- 
type junctions are present in the network. 





FIG. 5: Same as Fig. 3, except that now e — —0.4. In this case 
we start with initial conditions cho sen ra ndomly in the squar e 
whose vertices are the vacua (±^/15/7, 0) and (0, ±\/15/7). 
Here again only stable X-type junctions are present in the 
network. 



FIG. 7: The phase space distribution for the four timesteps 
of the simulation in Fig. 3. The vertices of the square in 
the plane (4>i,(j)2) are the minima of the potential, given by 
(±y/l5/7),±,/l5/7). 
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FIG. 6: Same as Fig. 3, except that now e = —0.8. Here we 
start with initial conditions chosen randomly in the square 
whose vertices are the vacua (±yl5/4, 0) and (0, ±^15/4). 
Note that in this case only K-type junctions are present in 
the network. 



FIG. 8: The phase space distribution for the four timesteps 
of the simulation in Fig. 4. The vertices of the square in 
the plane (<^>i,</>2) are the minima of the potential, given by 

(±\/l57l3),±yi57l3). 



only the edges are populated. These numerical results 
confirm our discussion above, as well as our original re- 
sults in Paper I. We emphasize that there is a valuable 
lesson to be learned at this point: the type of stable junc- 
tions that will form (in physical space) in a given model 
is an issue that is related to, but quite distinct from, 
the geometry of the various vacua in configuration space. 



In other words, there are both topological an energetic 
considerations involved in determining the detailed prop- 
erties of a given model, and one must be mindful of both 
of them in any rigorous analysis. 

Finally, we should also emphasize that the particular 
cases where er 2 = and er 2 = -1 do not represent 
marginal cases between the formation of the junctions 
of the type Y and X. Instead, these choices decouple 
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and writing the potential in terms of these. One then 
finds that these rotated fields decouple for er 2 = — 1. 

For an interesting example of a marginal case where 
both types of junctions are allowed, we can consider the 
model described by a complex scalar field $ with La- 
grangian 



C = d^d^ - V{§) 
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and the potential 



FIG. 9: The phase space distribution for the four timesteps 
of the simulation in Fig. 5. The vertices of the square in 
the plane (4>i,(j)2) are the minima of the potential, given by 
(±VT577).0) and (0,±vT5/7). 



V($) = K \<$> N - 1 ; 



(9) 



where k is a real parameter and N is integer. This model 
has vacua at 




<& = e L - 1 7i = 0, 1,...,7V-1. 

We can alternatively define the field's phase as <fi, and 
it is then obvious that we will have N minima, evenly 
spaced around </>. The case N — 2 produces standard 
domain walls and TV = 3 produces Y-type junctions, but 
the case N = 4 is slightly more subtle. The potential (9) 
has supersymmetric properties, and hence the energy of 
a specific solution depends only on the initial and final 
vacua. In other words, the possible ways of connecting 
two opposite vacua (directly or through the intermediate 
vacuum) will have the same energy. Hence this case is 
an example of the scenario where tr^ = 2er e . There is 
therefore no local energetic argument preferring one type 
of junction to the other, and consequently Y-type and 
X-type junctions will always co-exist. This can also be 
confirmed numerically, and visually the results can be 
approximately described as the superposition of the two 
distinct cases shown above. 



FIG. 10: The phase space distribution for the four timesteps 
of the simulation in Fig. 6. The vertices of the square in 
the plane (0i,(/>2) are the minima of the potential, given by 
(±^/l5/l),0) and (0,±y/TE/l). 



the two fields, and so no junctions are formed anymore. 
Consequently simulations will simply show the superpo- 
sition of two single fields evolving independently. This is 
trivial to see for the case er 2 = 0. For the case er 2 = — 1, 
this is best seen by defining a new pair of scalar fields ip\ 



III. THREE-FIELD MODELS 

We now move on to the case of models with three scalar 
fields. The discussion turns out to be quite similar to 
the two-field case, though at some points the effect of 
the increased dimensionality provides for different phe- 
nomenology. We will again consider the analogous BBL 
model [21], but we will also discuss the perturbed 0(3) 
Kubotani model, which has been studied in the past, 
starting with [5]. In fact we shall show that the two 
models are very similar, though not quite identical. 



The BBL model 



The extension of the BBL model of [21] for three real 
scalar fields is straightforwardly given by 



1 3 



V(M), 



(10) 



where the 
the form 



are real scalar fields and the potential has 
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where again r and e are two real parameters. As in the 
two-field case, there are two branches for the minima. In 
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form 



-2/5 < er 2 < 1/2 the minima are of the 
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while for — 1 < er 2 < —2/5 they are of the form 
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FIG. 11: Types of junctions (bottom) and configuration of 
the minima (top), as a function of the parameter er 2 in the 
two-field BBL model [21]. Notice that in the case er 2 > 
there are effectively two kinds of Y-type junctions, involving 
either the internal diagonal or the external diagonal sector. 



junctions which involve the axes sector (as opposed to 
only the edge sector), correspond to field space configu- 
rations where one of the three fields vanishes. The case 
e = is again not interesting because it represents the 
evolution of the three decoupled fields, with no junctions 
being formed. A numerical study can be performed along 
the same lines as what was done for the two-field case, 
and its outcome is summarized in Fig. 11. 
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In the former case there are 8 minima, which are placed 
at the vertices of a cube in the space {<j)\, <f>2,<f>3). I n the 
latter one, there are 6 minima, which are placed at the 
vertices of an octahedron. The minima of the second case 
can alternatively be thought of as being located at the 
centers of the faces of a cube. 

In the first case (12), there are twenty-eight topological 
sectors and three kinds of walls, which for obvious reasons 
we can refer to as edges, external diagonals and internal 
diagonals. The number of different walls of each type 
is respectively twelve, twelve and four. In the second 
case (13), there are fifteen topological sectors and two 
kinds of walls, which we can refer to as edges and axes. 
The number of different walls of each type is respectively 
twelve and three. 

Again the choice of the parameter e determines what 
type of junctions will be present. For the branch er 2 > 
—2/5, corresponding to the cubic solution of Eq. (12), 
X-type junctions survive only if the junctions involving 
walls from the diagonal sectors are energetically disfa- 
vored. This requires aid > 3er e (for the internal diago- 
nals) and a ec i > 2o~ e (for the external ones). These condi- 
tions are only verified in the interval — 2/5 < er 2 < 0, so 
in this range we do have stable X-type junctions. Con- 
versely, for < er 2 < 1/2 only the Y-type junctions sur- 
vive, and these may involve either of the diagonal sectors, 
so we effectively have two types of Y-junctions. On the 
other hand, for the branch er 2 < —2/5, corresponding 
to the octahedral solution of Eq. (13), both Y-type and 
X-type junctions can survive. Note that in this octa- 
hedron branch X-type junctions, as well as any Y-type 



B. The Kubotani model 

The perturbed O(N) models are a different class of 
models with N real scalar fields, which have been stud- 
ied in [5, 22]. Since the case that has been studied in 
more detail is N = 3, it is interesting to contrast it with 
the BBL model. In full generality, the model has the 
following potential 
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the model parameters should be such that £ + A > and 
£ + AA>0. 

A few particular cases are worth pointing out. Firstly, 
the case A = 0,£ > is the decoupled limit. This has 
minima <f> 2 = Q 2 . Again there are no real junctions: the 
different types of walls simply pass through each other 
without interacting, whereas walls of the same type in- 
tercommute normally. Still this is interesting as a test for 
numerical codes. Secondly, the case £ = 0, A > is the 
0(iV)-symmetric case. Here the minima are (j> 2 = n 2 /N. 
The N-th case has a total of 2 N minima, corresponding 
to the vertices of an N-dimcnsional hypercube. 

In the general case the minima are of the form 
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FIG. 12: The matter-era evolution of the domain wall network 
of a Kubotani-type model with a parameter choice A = 3/20, 
r\ 2 — 10/3 and £ = 1/12. This also mimics a three-field BBL 
model with r — \/3/2 and e = —0.2. Note that only stable X- 
type junctions survive. From left to right and top to bottom, 
the horizon is approximately 1/16, 1/8, 1/4 and 1/2 of the 
box size respectively. 



FIG. 13: The matter-era evolution of the domain wall network 
of a Kubotani-type model with a parameter choice A = 3/10, 
r/ 2 — 5/3 and £ = —1/6. This also mimics a three-field BBL 
model with r — \/3/2 and e = —0.4. Note that both Y-type 
and X-type junctions survive. From left to right and top to 
bottom, the horizon is approximately 1/16, 1/8, 1/4 and 1/2 
of the box size respectively. 



for the case £ < 0. The analogies with the BBL model 
should be fairly clear. In general the former has 2^ min- 
ima, whose vertices form an N-dimensional (hyper)cube, 
while the latter has 2N minima, that can be thought of as 
lying in the center of the faces of the said N-dimensional 
cube. The wall thickness is 
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while the wall tension is 
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Another particular subclass of models has £ = and 
N = 3; these were first studied in [5], and for this rea- 
son are usually referred to as Kubotani models, and we 
will use this choice for our numerical illustrations. The 
key difference between this class of models and the BBL 
models is the fact that for £ > the tension of the di- 
agonal walls (of which there are (N — 1) types in the 
iV-dimensional case) is always higher than N times the 
value of the edge walls. Consequently, for £ > the sta- 
ble junctions will always be of X-type. In other words, 
unlike the BBL model, Y-type junctions will never be 
stable in the cubic branch. The results of one simula- 
tion in this regime are shown in Fig. 12, where we have 
chosen A = 3/20, rj 2 = 10/3 and £ = 1/12. 

For £ < 0, it is again clear that Y-type junctions will 
form. But just like in the BBL case, one can also envisage 
having X-type junctions, which can arise from field space 
configurations where one of the three scalar fields van- 
ishes. Notice that in fact a chessboard- type lattice could 




FIG. 14: Two Y-type junctions in a Kubotani-type £ < 
model decay into a stable X-type junction; see the text for a 
detailed discussion. 



be built, with just four of the vacua lying on a plane, that 
might be stable. However, no realistic spontaneous sym- 
metric breaking mechanism can occur that selects only 
a subset of the vacua. Therefore the existence of both 
kinds of junctions is expected in any realistic simulation. 
This fact can be confirmed in Fig. 13 for the param- 
eter choices A = 3/10, rj 2 = 5/3 and £ = -1/6. We 
also note that even though the majority of the junctions 
are of Y-type, the fraction of X-type junctions remains 
approximately constant during the evolution. This issue 
deserves further study. 

This can be seen more easily if we construct by hand 
a simple square network containing four distinct vacua. 
In Fig. 14 we start with four vacua lying on a plane and 
forming two Y-type junctions: these are destroyed and a 




FIG. 15: An X-type junction in a Kubotani-type £ < model 
decays into two Y-type junctions; see the text for a detailed 
discussion. 



new X-type junction is formed. Conversely, in Fig. 15 
we have a set of four vacua that are not on a plane in 
field space: in this case any X-type junction is unstable, 
and decays into a pair of Y-type ones. 



C. Relating the BBL and Kubotani models 

The above discussion makes it clear that the BBL and 
Kubotani models are quite similar, but it is actually pos- 
sible to exhibit a direct relation between both potentials. 
The easiest way to see this is to relate them through the 
generalized potential 
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where A, B, C and D are real parameters. For the values 
for the parameters given by 
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one has exactly the potential of the BBL model. On the 
other hand, the choice 
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= A + e, B = -2\ V 2 , 


(22) 
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leads to the potential for the Kubotani model (recall that 
this has £ = 0). 

Using the relations above, it is easy to see that the 
models are related by 
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where A, £ and r\ are the Kubotani parameters and e, 
r, are the corresponding BBL parameters. The relations 
above lead us to conclude that for A > and ij 2 > 0, 
only the interval e < of the BBL model can be mapped 
into the Kubotani model. In addition, we see that £ > 




FIG. 16: The matter-era evolution of the domain wall network 
for a three-field BBL model with r = y3/2 and e = 0.1. Note 
that only stable Y-type junctions survive. From left to right 
and top to bottom, the horizon is approximately 1/16, 1/8, 
1/4 and 1/2 of the box size respectively. 



is related to the range —2/5 < er 2 < while £ < is 
related to —1/2 < er 2 < —2/5. However, notice that 
there are two consistency conditions, 2rj\ 2 = 1 on the 
Kubotani side, and er 2 = -27e 2 /2 - 2/3 on the BBL 
side (though the latter affects only the constant term in 
the potential), so even in this restricted e < range the 
correspondence is not one-to-one. 

Hence our Kubotani-model simulations in Figs. 12 and 
13 can also be interpreted as BBL model simulations. In 
Fig. 12 we took A = 3/20, r/ 2 = 10/3 and £ = 1/12. 
Note that only X-type junctions are formed. Using the 
relations above w e fin d that this corresponds to the BBL 
model for r = \/3/2 and e = —0.2, and similar results 
are obtained for the entire interval —2/5 < er 2 < 0. On 
the other hand, in Fig. 13 we have used A = 3/10, rj 2 = 
5/3 and £ = —1/6 in the Kubotani model (r = \/3/2 
and e = —0.4 in BBL model). Only Y-type junctions 
survive for any chosen negative value of £ (that is —1/2 < 
er 2 < —2/5). Finally, for completeness we present the 
case for < er 2 < 1/2 in the BBL model, which has no 
corres pond ence in the Kubotani case. We have chosen 
r = y/3/2 and e = 0.1 to illustrate this scenario in Fig. 
16. 



IV. THE IDEAL MODEL 

The availability of a large number of models (of which 
the two studied above are but examples), all different 
but to some extent related, begs the question of which 
features are fundamental and which are irrelevant de- 
tails. More to the point, one might ask what is the best 




FIG. 17: Illustration of a random domain distribution in a 
planar network taken from the fourth snapshot of Fig. 13. 
Each domain has x edges and each junction is formed by d 
walls. The average values are (d) = 53/17 and (x) = 106/19. 



possible model for domain walls, at least from the point 
of view of its potential to produce frustrated networks. 
Based on our results so far, we shall discuss what the key 
characteristics of such a model are, and discuss an ex- 
plicit construction. It turns out that a model considered 
in [23] has some of the characteristics of this ideal model. 
In order to make the argument clear, we shall go through 
the key steps in considerable detail, starting with a dis- 
cussion we already presented in Paper I and building up 
from there. 

In a two-dimensional domain wall network generated 
dynamically by a cosmological phase transition, the num- 
ber of edges of the domains can assume several values, 
with a particular distribution. The dimensionality of the 
junctions in such a network, that is, the number of walls 
meeting each junction, depends on the structure of the 
chosen potential. One example is the Kubotani model 
that can be seen in Figs. 12(for £ > 0) and 13 (for 
£ < 0). While in the former case only X-junctions are 
present, the latter one presents both Y-type and X-type 
junctions. Note in both cases the stable junctions are 
formed by walls with equal tension. We should empha- 
size that as the average dimensionality of the junctions 
increases, the average number of faces of the domains will 
be smaller. 

As a simple illustration, we can take the fourth snap- 
shot in Fig. 13, which is schematically shown in Fig. 17. 
The average number of faces x of each domain is related 
to the dimensionality d of the junctions by 



«a:>-2)«d>-2) =4. 



(25) 



In Fig. 17, these values are (d) = 53/17 ~ 3.11 and 
(x) = 106/19 ~ 5.58. We have also verified that these 
numbers are approximately constant during the evolution 
of the network. Indeed, one can in principle use these 
numbers to characterize each particular model, that is 
each choice of the potential and of the corresponding free 
parameters. This is another issue that is beyond the 
scope of the present discussion, but to which we shall 
return in the future. 



Now, the stability of each particular domain will de- 
pend, among other things, on its number of edges. Let 
us momentarily assume that d is fixed to be equal to 3. 
By using the Euler relation in Eq. 25, such network must 
have (x) — 6. Of course the network will contain domains 
with different values of x. However we have shown in Pa- 
per I that, in this case, only domains with x > 6 might 
possibly be stable, and that all domains with x < 6 will 
quickly collapse in order to minimize the energy. There- 
fore, no domain wall evolution starting by random initial 
conditions in such models can provide a stable network. 
Furthermore, in scenarios where there are domain walls 
with different tensions, an additional source of instabil- 
ity is also present: when these heavier walls collapse, the 
dimensionality of the corresponding junction is increased. 

Putting all these facts together, the logical conse- 
quences are clear. In an ideal model, the probability that 
two domains in the same vacuum state are close to each 
other should be as small as possible. This can be accom- 
plished in models with a very large number of different 
vacua (say, having the number of scalar fields involved 
N — v oo). In this case the collapse of a single domain 
will only very rarely lead to the fusion of two of the sur- 
rounding domains. Note that this fusion would lead to a 
further reduction of the total number of edges of the con- 
tiguous domains which as we saw is clearly undesirable 
since it would increase the probability that some of the 
contiguous domains would themselves become unstable 
to collapse. On the other hand, all the possible domain 
walls should have equal tensions. Again, if that were 
not the case we would be adding a different source of in- 
stability since the walls with higher tension would tend 
to collapse, thereby increasing the dimensionality of the 
junctions which, in turn, would lead to a decrease in the 
average number of edges of individual domains and to 
the production of further unstable two-edge domains. 

The ideal model (as far as frustration is concerned) is 
therefore a model with a very large number of vacua with 
all the domain walls connecting the various vacua hav- 
ing the same tension. In Paper I we have already shown 
that, due to energetic considerations, the stable junctions 
of such a model must necessarily be of Y-type only, so 
that d = 3. We have also shown that the only possible 
equilibrium configuration in two dimensions with Y-type 
junctions is the 'honeycomb' hexagonal-type lattice (oth- 
erwise unstable two, three, four and five edge domains 
would necessarily occur). The fact such lattices are not 
really stable (they allow for locally confined energy defor- 
mations although the Hubble damping may temporarily 
prevent the collapse of such configurations) together with 
the very small probably that such a lattice could be an 
attractor for the evolution of domain wall network simu- 
lations from realistic initial conditions has led us to our 
no frustration conjecture. However, the analysis of Pa- 
per I can only be fully applied to domain wall networks 
in two spatial dimensions, and therefore the construction 
of an explicit realization of the ideal model is crucial if we 
want to test our conjecture in the more realistic case of 
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FIG. 18: The illustration of the vacua displacement of three 
Higgs potential. Note that there are only three equidistant 
points where these minima coincide. 




FIG. 19: Possible collapses of a domain with four edges in the 
model with only three vacua. 



three spacetime dimensions. It is this missing step that 
we now provide. 

Geometrically speaking, the ideal potential is there- 
fore described by N real scalar fields and has mutually 
equidistant minima. The number of vacua is AT + 1 for 
N real scalar fields, and the energetic cost for a specific 
transition between any two of them is the same. An ex- 
plicit realization of such a model is given by the potential 



V = 



X 



N+l 



N- 



iE'SW- 



with 



N 



= £to-^) 2 



(26) 



(27) 



where pi . are the N + l coordinates of the vacua of the 
potential. We have chosen p$, to be the vertices of an 
(N + l)-dimensional regular polygon, and fixed the dis- 
tance between the vacua to be equal to the parameter vq. 
The model is therefore the sum of N + 1 <fi 6 potentials. 
Each of them has one minimum located at the center, 
and a continuum set of minima at a distance tq from the 
center. Note that N of these vacua are located exactly 
at the centers of the other potentials. 

For instance, if one takes N — 2 here we will have 
three equidistant vacua at the vertices of an equilateral 
triangle, as illustrated in Fig. 18. This particular case is 




FIG. 20: Possible collapse of a domain with four edges in the 
model with more than five vacua. 



<5<3 



FIG. 21: The vacua of the Carter model [23] 
two phase angles <j>\ and <f>2. 



in terms of the 



in fact analogous to the model given by Eq. (9) if one 
takes N = 3 in that potential. Despite the guarantee 
that the only stable junction will be of Y-type, we do 
not expect that a stable lattice will emerge from random 
initial conditions. Since this model has only three vacua, 
the collapse of domains with four edges will always hap- 
pen as illustrated in Fig. 19. In some sense, this N = 2 
realization of Eq. (26) is maximally unstable. It is obvi- 
ous that if this potential had two more minima, the decay 
shown in Fig. 20 would become possible. In this config- 
uration the collapse of the square occurs decreasing the 
number of edges of the neighbor domains, but not cutting 
the network. Notice that the probability of appearance 
of such a configuration will still be very small, though en- 
hanceable by the addition of further vacua. Again, this 
is the reason why we want to have N — ► oo. 

The case with A~ = 3 has four vacua located at the 
vertices of a tetrahedron. All the minima are connected 
to each other by walls with equal tension. We emphasize 
that although the two-field BBL model also has four min- 
ima, the two models are completely different. In the BBL 
case stable X-type junctions may be formed depending 
on the parameter e, and when there are Y-type junctions 
two different kinds of wall exist, with different tensions. 
Here only Y-type junctions are formed and all of them 
have the same tension. 

Let us point out a final example, A~ = 4, which has four 
fields whose five minima are located at the vertices of a 
pentahedron. A very similar analogue was introduced by 
Carter in Ref. [23]. This indeed has five evenly spaced 
minima, and all the walls are of the same energy — the 
only difference is that Carter's model is of <jr type, while 
ours is of </> 6 type. Obviously, both this model and ours 
allow for the configuration in Fig. 20. However, as we 
have discussed, this configuration is not enough for en- 
suring the stability of network frustration. Note that we 
again see that Carter's claim that this model will form 
stable X-type junctions is manifestly incorrect — we will 
also confirm this numerically below. 
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Note that it is also possible to provide a generalization 
of the Carter model, for the case of a generic (though 
even) number of real scalar fields. In [23] it is pointed 
out that the model can be thought of in terms of ei- 
ther four real scalar fields or two complex scalar fields. 
Indeed, the phases (let's call them ipi and ^2) of these 
two complex scalar fields provide the simplest way of vi- 
sualizing the configuration of the vacua — the alternative 
of (N + 1) evenly spaced points on and (N — l)-sphere 
being obviously not very clear. Since the minima must 
be evenly spaced, they will lie on the surface of some 
sphere, and one can therefore move between them by 
moving along some suitably chosen phases. For N = 2M 
real scalar fields (M complex fields) there will therefore 
be M such phases. Let us consider each such phase indi- 
vidually: it will have a certain number of vacua, equally 
spaced along it, but obviously only two nearest neighbors. 
But since there are M such phases, we immediately see 
that the total number of equidistant minima should be 
2M + 1 = N + 1. In Fig. 21 we can see the configu- 
ration of the minima for the Carter model [23] . Notice 
that starting at any of the minima, we can reach any of 
the other minima in a single direct displacement, either 
moving along -0i or along fa as the case may be. 

A generalization of the potential given in [23] , which is 
always of 4 type, is the perturbed [£7(1)] M model which 
has the following form 



M 



M 



2/M 



v(*0 = £(*? 



-£ 2 /m( 



II*? 



(28) 



Here 3>j are the M complex scalar fields, with </>, the 
corresponding phases, and /m(0z) 1S a function of these 
phases, having 1M. + 1 minima and taking the value 
fMmin = —M. at those points. The minima of the poten- 
tial will then be of the form 



& 



>r 



l-e 



(29) 



and be evenly spaced for a suitable choice of the function 
Jm- Specifically, for M = 1 we trivially have 



/i = cos(30i), 



(30) 



yielding the three evenly spaced minima we already men- 
tioned, while for M = 2we should choose 



f 2 = cos (2^1 + fa) + cos (202 - 0i ) , 



(31) 



which is the potential discussed in [23] , and yields the five 
evenly spaced minima shown in Fig. 21. The construc- 
tion of generalizations for higher values of M is straight- 
forward in principle, though algebraically tedious in prac- 
tice. 

To confirm our description of the model (26) we have 
simulated a number of realizations of the ideal class of 
models, that is the potential Eq. (26) for different values 
of N up to N = 10. In all these simulations we have 




FIG. 22: The matter-era evolution of the domain wall network 
for the ideal model with N — 4. Note that this model is 
effectively analogous to that of [23], and that only stable Y- 
type junctions survive. From left to right and top to bottom, 
the horizon is approximately 1/16, 1/8, 1/4 and 1/2 of the 
box size respectively. 



confirmed that only stable Y-type junctions are formed. 
In addition, any X-type junction imposed in the network 
by hand have been destroyed quickly. Two examples are 
shown in Figs. 22 and 23, respectively for N = 4 and 
N = 7. Note in particular that these results confirm that 
in the Carter (pentahedral) model only Y-type junctions 
survive. It is very interesting to study several network 
properties, such as the degree of instability of the net- 
work (along the lines outlined above) or the degree of 
convergence towards a scaling solution [16-18] as a func- 
tion of the number of fields N. This is again a task that 
is left for future study. 



V. CONCLUSIONS 

We have studied several typical domain wall forming 
models in two spatial dimensions, and discussed the con- 
ditions under which various types of defect junctions can 
exist. In particular, we focused on the BBL [21] and Kub- 
otani [5] models, both of which have been studied in the 
past. We have highlighted the situations where Y-type 
junctions are stable, those where only X-type ones are 
stable, and those where both types co-exist. Using both 
analytic and numerical techniques we have shown that 
topological and energetic arguments both contribute to 
the types and behavior of the allowed junctions, and we 
have identified their respective roles. A through under- 
standing of the effects of both of these mechanisms is vital 
for a proper assessment of the cosmological consequences 
of defect models with junctions. This is particularly true 



12 




FIG. 23: The matter-era evolution of the domain wall network 
for the ideal model with N = 7. Note that only stable Y-type 
junctions survive. From left to right and top to bottom, the 
horizon is approximately 1/16, 1/8, 1/4 and 1/2 of the box 
size respectively. 



for the case of cosmic (super)strings, which haven't as 
yet been explored in quantitative detail. 

Our analysis also allowed us to distinguish which fea- 
tures of a particular model are crucial for its behavior 
and which ones are, to some extent at least, irrelevant 
details. The crucial step in the analysis is the simple but 
previously neglected result (which we have discussed in 
Paper I) that in models with more than two vacua, if all 
the domain walls connecting the various vacua have equal 
energies then only Y-type stable junctions can be stable. 
This in turn allowed us to define and explicitly construct 
a realization of what would be the ideal model in terms 
of would-be frustrated networks. Granted that this is a 
toy model, it does have the advantage of including all 
of the features that are desirable in terms of frustration 
(and none of the undesirable ones), as well as the further 
advantage of being simple enough to be easy to explore 
numerically. This model is of 4> e type, and our construc- 
tion will work for any number of real scalar fields, as 
well as for any number of spatial dimensions. In addi- 
tion, we have also provided a different realization, which 



is to some extent simpler (being of 4> 4 type) but arguably 
less general, since it applies only to an even number of 
scalar fields — though one could of course claim that it 
applies to any number of complex fields. 

These results strongly support our no frustration con- 
jecture. In all our 2D numerical simulations, including 
those of the ideal model, we find a dynamical behavior 
consistent with a convergence towards a scaling solution 
[16-18] and no evidence of any behavior that could even- 
tually lead to frustration. (We also note that recent work 
by Battye et al. provides further support for this.) We 
therefore conclude that in two spatial dimensions no frus- 
trated network can be formed dynamically out of random 
initial conditions that would mimic a realistic phase tran- 
sition. Notice that we are making no claims about the 
possibility of designing (by hand) a lattice-type configu- 
ration that could (in rather particular and/or idealized 
cases, such as the static limit) be stable against small 
perturbations. That is of course an interesting mathe- 
matical problem, but whether or not such a possibility 
exists is of very little consequence for the cosmological 
behavior of this type of defects. 

What now remains to be done is to extend this anal- 
ysis to the more realistic case of three spatial dimen- 
sions. With our identification of the ideal class of models, 
this means are now available to undertake this analysis. 
Notice that the increased dimensionality will have non- 
trivial effects in the behavior of the junctions. Neverthe- 
less, we have reasons to believe that the no frustration 
conjecture will also hold in three spatial dimensions. Wc 
leave a detailed analysis of this issue, as well as that sev- 
eral other interesting questions that we have highlighted 
at various points of this paper, for forthcoming work. 
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